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Abstract 

A superconducting ring has different sectors of states corresponding to differ- 
ent values of the trapped magnetic flux; this multitude of states can be used 
for quantum information storage. If a current supporting a nonzero flux is set 
up in the ring, fluctuations of electromagnetic field will be able to "detect" 
that current and thus cause a loss of quantum coherence. We estimate the 
decoherence exponent for a ring of a round type-II wire and find that it con- 
tains a macroscopic suppression factor (5/Ri) 2 , where R± is the radius of the 
wire, and 5 is the London penetration depth. We present some encouraging 
numerical estimates based on this result. 

PACS numbers: 03.67.Lx 
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I. INTRODUCTION AND RESULTS 



A superconducting ring is a good candidate for macroscopic quantum memory [l]. Its 
ability to store quantum information stems from the complicated "vacuum structure" , viz. 
the presence of different sectors of states corresponding to different values of magnetic flux 
through the ring. 

It is easy to suppress spontaneous transitions between these different sectors by making 
the ring sufficiently thick M. In the absence of such transitions, the main mechanism by 
which quantum information stored in the ring deteriorates is decoherence, i.e. a decrease, in 
time, of matrix elements between states with different values of the flux. Decoherence is due 
to interaction of the macroscopic variable, in our case the supercurrent in the ring, with the 
"environment" comprised by local degrees of freedom — fluctuations of the electromagnetic 
field. 

The supercurrent density required to support a single flux quantum through a macro- 
scopic ring is very small. This makes that state locally almost indistinguishable from the 
state with no flux at all. Based on this observation, we have suggested |J that in this case 
decoherence may not be as strong as it is generally believed to be for macroscopic systems. 
The purpose of the present paper is to discuss decoherence in the ring quantitatively. 

In general, the environment starts out in an arbitrary initial state with some density 
matrix pq. Two different values of the macroscopic variable give rise to two different histories 
of the environment, which are represented by evolution operators Ui(t) and U^if). We define 
the amount of coherence (at time t) between these two histories as the "overlap" 

C(t)=Tr[U 1 (t)p Ul(t)]. (1) 

Next, we define a decoherence exponent D(t) by 

exp[-D(t)} = \C(t)\ 2 . (2) 

For an environment comprised by linear oscillators, which in addition interact linearly 
with the macroscopic variable, decoherence had been widely discussed in connection with 
macroscopic tunneling |2],|3[]. That problem is similar but not identical to ours: in quantum 
memory tunneling is suppressed. Indeed, our problem is more analogous to the "error-free" 
case of ref. j|; however, we will be able to supply definite values for all of our parameters. 

In principle, the linear approximation can break down at large times even for weakly 
nonlinear environments. In addition, whatever we designate as the environment will inter- 
act with a still larger system and will loose energy to that system at some characteristic 
dissipation rate. For example, in our case, the electromagnetic field can be absorbed on the 
surface of the superconductor, as well as on the surfaces of any surrounding objects. In such 
cases, the evolution of the environment can no longer be represented by unitary operators. 

It is certainly of interest to understand how nonlinearity and dissipation in the environ- 
ment affect quantum coherence in the original system. However, for our present purposes 
the linear approximation will be sufficient. A mechanism for decoherence exists already at 
the linear level Q . Dissipation will only be discussed phenomenologically, to make sure that 
it does not significantly alter the evolution of the oscillators at relevant times. 

Specifically, we consider the following problem. Suppose at the initial time t = a su- 
percurrent has been set up in the ring but the electromagnetic field is still in the thermal 
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state corresponding to zero current, with a Gibbs distribution w m of probabilities for the 
energy eigenstates \m). Because the field is displaced from its true equilibrium (correspond- 
ing to a nonzero average magnetic field), it will oscillate. Consider two oscillation histories, 
corresponding to two different values of the current (one of which may be zero), and call 
their evolution operators Ui(t) and U^if)- The decoherence exponent in this case will be 
determined by 



exp[-D(t)]=J2\HU j 2 U 1 \m) w m . (3) 

m 

We describe a calculation of D(t) below. 

The problem just posed is admittedly somewhat artificial. Imagine that we switch from 
the state with no flux (and no current) to the state with unit flux by dragging a single flux 
line through the bulk of the superconductor. There is no reason to expect that the current 
will equilibrate first, and the field next, as in our model problem. Nevertheless, we think 
that it is good to know the answer to this model problem before turning to more realistic 
setups. 

The answer is this. Consider a circular ring of radius Rq made out of a London-type 
superconducting round wire of radius R±, in the limit 

5 < Ri < R , (4) 

where 5 is the London penetration depth. We find that at intermediate times R\/c <C t <C 
Rq/c (c is the speed of light), the exponent D(t) grows as t for t <C h/k B T and as t 2 for 
t ^> h/k B T, up to some logarithms. At large times, t 3> Ro/c, the exponent saturates at a 
certain limiting value, for which we obtain the following estimate: 

d f{u) ( nc \ (-Y W 

hm 16a EM ln 4 (i? /i?i) U o; mi J [ft J ' 1 1 

where «em is the fine-structure constant, and u = k B T '/hu m i n with u m i n ~ c/Rq; f{u) = 1 
for small u, and f(u) = u for large u. 

The main result is the macroscopic suppression factor (S/Ri) 2 in ([5]). Taking for estimates 
5 ~ 10~ 5 cm, Rq — 1 cm, R 1 = 1 mm, u min = nc/R Q ~ 10 11 s -1 , and replacing the 
logarithm (which controls logarithmic accuracy) with unity, we obtain x ~ T/(l K) and 

Ainx-lO" 7 /^)- 

These results apply not only to a ring made out of a solid superconductor, but also 
to a ring coated with a superconducting film, as long as the thickness of the film is much 
larger than the London depth 5. Similar calculations can be done for any type of persistent 
circular current, of which a supercurrent in a macroscopic ring is but one limiting case. The 
opposite limit is a circular current in a single atom — e.g. a circular Rydberg state ||, which 
has already been proposed as a basis for quantum computation . It would be interesting 
to do the calculation for that case and see how the decoherence rate changes as the current 
loop becomes microscopic. 

We now outline the main steps leading to estimate @. 
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II. THE OSCILLATOR HAMILTONIAN 



We will only consider decoherence at times 

t > Ri/c , (6) 

when, as we will see, the main contribution to the decoherence exponent comes from elec- 
tromagnetic modes with low frequencies: 

u < c/Ri . (7) 

For such modes, the response of a superconductor of the London type (the only type we 
consider) is determined, to a good accuracy, by the London penetration depth S. 

To count the oscillator modes, we enclose the ring in a large normalization sphere of radius 
1Z 3> Ro, centered at the center of the ring. We define a system of spherical coordinates 
(r, 9, 0), with angle 9 measured form the axis of the ring's rotational symmetry. 

Only the azimuthal component of the vector potential, indeed only the ^-independent 
modes of it, interact with the current. For these modes, the vector potential takes the form 

A{r,6,<f>]t) = e 4 ,{<f>)A(r,e ] t), (8) 

where is the azimuthal unit vector. Then, 

V 2 A = e, (V - — U-) A , (9) 
V rsin 9 J 



and the wave equation for A takes the form 

A - VM + ' 



r 2 sin 2 9 



V{r,9) 



A=^J. (10) 



Here j is the constant (in time) supercurrent density, and V(r, 9) is the "potential" that 
represents the Meissner effect: inside the superconductor 

V(r,0) = l/6 2 , (11) 

while outside the superconductor both j and V are zero. 

The equation of motion ( 1Q|) can be obtained, after substitution (H), from the effective 
Hamiltonian 



H = KM^ A + + ^ A 2 " 7 A j ) > (I 2 ) 



2 1^2 V a2 1 

l + — B + — A - -. 

■57rc z 8n 8n c 

where B = V x A is the magnetic field, and j = e^j. Such Hamiltonian formulation neglects 
any dissipation (absorption) of the electromagnetic field. This approximation is justified in 
Appendix A. We stress again that decoherence we calculate in this paper is not due to 
any such dissipation, but is a result of a "measurement" of the supercurrent, performed by 
fluctuations of the electromagnetic field. 

Now, expand the field A in a complete set of modes f n (r, 9): 
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A(r,9;t) = J2Ur,9)X n (t) . (13) 

n 

This expansion defines a set of oscillator coordinates X n . The modes /„ satisfy the following 
equation: 



■V 2 + ^^ 7j +V(r,6) 
r z sin t) 



J 1 

fn n~ fn j (1^) 



where uj n are real eigenfrequencies, together with the continuity conditions at the surface 
of the superconductor and some boundary condition at the boundary of the normalization 
volume. They are normalized in that volume by 

d?xf n (r,6)f n ,{r,6) = <W • (15) 

A typical low- frequency mode is shown schematically in Fig. 1. Notice that its main support 
is outside the superconductor. 

Substituting the expansion (|13"D in the Hamiltonian (|T2| ) we obtain the Hamiltonian for 
the oscillators: 

#=^E(*n + U 2 n Xt) "-E X « / ^3 (r, )fn(r, 0) , (16) 
oTCC n C n J 

which is of the Caldeira-Leggett || type. The role of a macroscopic variable is played by 
the total supercurrent I. Because in our case the total supercurrent has no dynamics of its 
own, it is not necessary to split each integral in (]l£j) into a product of q = / and a coupling 
C n . Nevertheless, to conform to the established notation, we will denote these integrals as 

qC n = — fd 3 xj(r,d)f n (r,9). (17) 



The quantum-mechanical problem specified by the Hamiltonian ([T6|) is exactly solvable, 
and the overlap (|l|) is easily calculable. For the decoherence exponent of two states corre- 
sponding to two different values q = qi and q = q2, with the initial state being a thermal 
ensemble at q = 0, we obtain 

D(t) = [°° duo^-\l - coscut] coth(ft/W2) , (18) 

nn Jo uj 

where Aq = q% — q 2 , (3 = 1/ksT, and J(u) is the spectral density defined as in refs. @||: 

m = ?E ^r 6 ^ - ») ; ( 19 ) 

according to (0) the masses of all oscillators are 

= *s? ■ (20) 

Eq. ([18D is similar to expressions found in the literature 0, but we should stress that ours 
is not a macroscopic tunneling problem: as we pointed out before, in quantum memory, 
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we strive to prevent spontaneous transitions (tunneling or otherwise) between states with 
different values of q. These different values now simply label various macroscopically distinct 
histories. 

Expression (|i~8D applies at both zero and finite temperatures. Two limiting cases can be 
considered: ksT <C huj, and fc#T ^> huj, where uj is a typical frequency at which the integral 
in (p~8|) saturates. In the first case ("low" temperatures), 

D{t) ^JAf_ fd^l-coBut]; (21) 
nh Jo uj 

in the second case ("high" temperatures), 

D ( t ) = r^^[l -cosut] . (22) 

nh z /3 Jo u 3 

To complete the calculation of D we now need to find the couplings C n , for the low- frequency 
modes of interest. 



III. DECOHERENCE AT INTERMEDIATE TIMES 

Our strategy will be as follows. We first compute the spectral density J{ui) for uj in the 
range 

c/R < uj < c/Ri . (23) 

The form of the spectral density at these uj determines the decoherence exponent at inter- 
mediate times 

Ri/c < t < R /c . (24) 

Then, in the next section, we establish that at smaller u, u c/Rq, J{uj) is proportional 
to uj 3 . According to (|i~8l) , that means that at t ^> Rq/c the exponent D(t) does not grow 
beyond a certain limiting value D\ im . Finally, we estimate -D nm . 
For each mode of oscillation, we define a wavenumber k n by 

k n = uj n /c. (25) 

At any time in the range (|2~4]), the main contribution to the decoherence exponent comes 
from modes with 

kn~l/ct. (26) 

For counting these modes, the normalization volume can be replaced by a tube, coaxial with 
the wire, of radius -Rnorm that satisfies 

Ct < .Rnorm < R ■ (27) 

This new normalization volume is more convenient than the original sphere because the 
problem now acquires an approximate cylindrical symmetry: we can view the wire as a 
straight cylinder, with fields subject to periodic boundary conditions at the ends. 
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This cylindrical problem has a natural set of cylindrical coordinates, (p,tp,z). The ex- 
pression (g) for the vector potential can now be approximated by 

A(p,<p;t) = e z A{p,ip;t) , (28) 

and the mode equation ( |14|) by 

-V 2 + V(p)]f n (p,p) = ^f n (p,p). (29) 



As seen from (|29"|), the mode functions f n are eigenfunctions of —id /dip, the angular mo- 
mentum about the z axis. We concentrate on functions with zero angular momentum (i.e. 
those independent of the angle <p), because for other values the coefficients C n in ( |19"D are 
suppressed by additional powers of k n R 1 . Thus, we are left with f n that depend only on the 
radial coordinate p, the distance to the wire's axis. 

By assumption, the wire's radius Ri is much larger than the penetration depth 5. Then, 
the solution to equation (fTl) inside the superconductor is approximately 

f n (p)^A n exp[(p~R 1 )/5] . (30) 

On the outside, the solution is a combination of Bessel functions; near the boundary, for k n 
satisfying (PB[), it reduces to 

f n (p) « A n + B n hx{p/R l ) . (31) 

Eq. ( pip already takes into account the continuity of the solution on the surface; the 
continuity of the first derivative requires B n = A n Ri/6, so (|31|) becomes 

/^)«A*(l + ^ln(p/#i)) . (32) 

Notice the enhancement of the coefficient of the logarithm; it comes from matching the first 
derivatives. 

The coefficient A n in (0), (|32|) is determined from the normalization condition (|l"5|) . To 
compute the normalization integral, we need the solution further away from the surface of 



the superconductor, where the form (|32|) no longer applies. But (|32|) does allow us to pick up 
the coefficients of the Bessel functions. Writing ln(p/i?i) = ln(£; n p) — ln(/c n i?i), we observe 
that the coefficient of Yo(k n p) (the function of the second kind) is proportional to A n R\/6, 
while the coefficient of Jo(k n p) (the function of the first kind) is additionally enhanced by 
ln(A; n i2i). 

We will be content with logarithmic accuracy; to such accuracy, away from the surface 
we have 

R 

fn{p) = -An-^-l^knR^Joiknp) . (33) 

To count the modes, we use the Dirichlet boundary condition at p = R n0 rm'- 

fn(Rnorm) = . (34) 
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Normalizing ([33|) in a cylinder of radius -Rnorm and length L, we find 



An -i ^ y /2 8 1 (35) 

V 2LR novm I R\ I In k n R\ 



Finally, using for j the equilibrium current density 



j(p) = 2^ eXp[(p_jRl)/5] ' (36) 



where / is the total supercurrent, we obtain for the spectral density 

(37) 

We now specialize to the case when q\ = 0, and q2 = I s is the current corresponding to 
one quantum of flux through the ring: to logarithmic accuracy 

C$n (38) 



2L\u(L/Rx) ' 
where $o = nhc/e is the flux quantum. We obtain 

,. , 2 T , . 7r 3 c 2 ^ 2 / 5 \ 2 1 . . 

{Aq) J{U) = l6e 2 Lln\L/ Rl ) [K ) lA^Rjcj (39) 

(e is the absolute value of the electron charge). 

With this form of J (a;), the integral in eq. (|18D gets the main contribution from u ~ 1/1 
So, at t <C h/k B T, we apply eq. fl21j), to obtain (to logarithmic accuracy) 

W ~ 16e 2 L In 2 (L/RJ \~R~J In 2 (ct/^x) ' ( ^ 

At i ^> h/kgT, we apply the "high" -temperature limit (1221). Although the integral in 



( p2|) in our case is convergent, we impose a low-frequency cutoff at = c/ i? ; a representative 
value below which eq. (|37|) is not applicable. This will give us an estimate of the accuracy 
of the result. We obtain 



D(t) 



t 2 n 2 c 2 k B T (jT 
8e 2 L\n 2 (L/R 1 ) \R lf 



2 r 



1 + o 



]n(ct/Rx) \n(R /R 1 ) {bi^ct/R^ 



(41) 



The second term in the bracket is a correction due to the low-frequency cutoff. Under the 
condition (fj4|), this term is small compared to the first term (the one from uj ~ 1/t), but as t 
approaches values of order Rq/cH becomes comparable to the first term. Thus, at t ~ Rq/c 
the expansion indicated in (|4l| ) breaks down, and we need a different method to establish 
the form of D(t) at large t. 
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IV. DECOHERENCE AT LARGE TIMES 



We now turn to modes with 

K < l/Ro • (42) 

For counting these, we use the original normalization sphere and the Dirichlet boundary 
condition 

f n (K) = . (43) 

The modes are normalized in the sphere by the condition (I5p. 

Let us expand each mode function f n in associated Legendre functions P t L (cos9): 

f n (r,9) = (2nN l )- 1 / 2 Y^S nl (r)P^cose) , (44) 

i>i 

where S n i are some radial functions, and Ni is the normalization integral for 

n, = 2 4±v (46) 

21 + 1 K J 

Note the absence of I = 0: there is no (nonzero) P . 

Eq. ([i~5|) now translates into the following normalization condition for the radial functions 

Sni' 

V / S n i(r)S nn (r)r 2 dr = b nn < , (46) 
i>i Jo 

while the mode equation (|14]) leads to the following equation for S n i: 

- ~ [r 2 ^Sn^ + l ^~^S nl - klS nl + ^2 V w S nV = . (47) 

The "potential" Vu> here is given by 

V lv (r) = (N^,)- 1 / 2 r V(r, 9) PH cos 9) Ph( cos 9) sin9d9 , (48) 
Jo 

in terms of the "potential" V(r, 9). 

The potential ( ^8l) is nonzero only for R in < r < R OVLt , where i? in and R out are the inner 
and outer radii of the ring. In the region r < R in , eq. ( f4"T| ) is a free radial equation, in which 
we can now neglect k\. A solution regular at r = is then 

S nl (r) = D nl r l , (49) 

where D n i are some constants. In the outside region, r > -R ou t, where eq. (f47D is also free, 
the solution is a linear combination of two spherical waves. As long as r l/k n , so that k\ 
is still negligible, we can write 
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S n i(r) = F nl ( -^-V + G nl ( ^V +1 , (50) 



R 



out 



with some constants F n i and G n [. On the other hand, at large distances r ^> 1 /k n , we have 
(cf. ref. 0) 



S ^ ~ + sinfor - ttZ/2 + fr.) , (51) 

rk l + L R l ont cosd n 



where 5 n is a phase shift. 

Because the smallest possible I is I — 1, from the normalization condition fljE| ) we see 
that 



O(^) (52) 



(the phase shift 5 n is O(^)). For each n, the sets of coefficients F n i, G n i, and D n i are related 
through scattering on the potential Vu>. Because this scattering occurs at r ~ _R <S l/fc n , 
it is insensitive to so G n i and _D n z, and therefore also /„ on the surface of the wire, are 
all O(kl). Then, according to eq. (|T§D, the spectral density J{oj) is 0(cu 3 ). 

This suppression of the spectral density at small u can be explained by noting that the 
large wavelength fluctuations of the field "see" simultaneously two diametrically opposite 
segments of the wire, with currents that add up to zero. Such fluctuations therefore have 
only derivative interactions with the current; hence the extra powers of k n , or uj. 

From eq. (|TED, we now see that the low frequency region u -C c/ Rq does not significantly 
contribute to decoherence, and the main contribution at t ^> Ro/c comes from modes 
with uj ~ c/Ro. We can estimate decoherence at these large times by using (|I^) with 
the expression (p9|), which correctly describes higher frequencies, and an infrared cutoff at 
^min ~ c/Rq. In this way, we arrive at the estimate (||). 

APPENDIX A: THE ROLE OF DISSIPATION 

Here we list some estimates for absorption rates of low-frequency electromagnetic field on 
the surface of a superconductor. The purpose is to show that the corresponding dissipation 
time is much larger than the timescale Ro/c, at which the decoherence exponent saturates. 

The response kernel Q, which determines the current induced in the superconductor by 
the electromagnetic field, is defined by 

Jind (x, t) = -J <Px!dt'Q{x, x'; t - t')A(x', . (Al) 

The low-frequency expansion of the Fourier transform of Q (inside the superconductor) is 

QH = ^-w/c 2 + ... , (A2) 



where the first term is the one taken into account in eq. (jT2|), while the second term describes 
dissipation due to a finite conductivity a. 
We only consider low frequencies, 
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uj <C c/Ri , 



(A3) 



and conventional (s-wave) superconductors with critical temperatures T c of order of a few 
Kelvin. For Ri ~ 1 mm, we have Tiuo <C 2A , where A is the zero-temperature gap. 
So, at sufficiently low temperatures, a single photon cannot break a Cooper pair, and the 
dissipation is entirely due to thermally excited quasiparticles. Therefore, we estimate the 
conductivity as 

a ~ a N exp[-A(T)/k B T] , (A4) 

where is the normal-state conductivity at T around T c , and A(T) is the temperature- 
dependent gap. 

Using ( |A4| ) with aN ~ 10 18 s _1 (corresponding to 10 8 f2 _1 m _1 in SI) and 5 ~ 10~ 5 cm, 
we find that for our range of frequencies the second term in ( |A2|) is much smaller than the 
first. This allows us to expand the surface impedance as 

C(w) = u[-A7tc 2 Q(uj)]- 1/2 = -i— [l + 27r^a5 2 /c 2 + 0(^ 2 )1 . (A5) 
The dissipative effect (absorption) is represented by the real part of £, 

( r {uj) = 2tiuj 2 (x5 3 /c 3 ~ 10- 6 (cj/10 n s _1 ) 2 exp[-A(r)/A;BTI » ( A6 ) 

where the estimate is for the same values of aN and S as before. 

For example, suppose that the device is inside a cavity, and the cavity and the ring are 
made of the same superconducting material. Then, assuming that radiation does not leak 
out, dissipation is mainly due to absorption on the walls of the cavity. The rate of dissipation 
is estimated as 

T- 1 ~ CW^cav , (A7) 

where i? cav > Ro is the radius of the cavity. The timescale r is thus 1/Cb times larger than 
-Rcav/c, which in turn is larger than Rq/c, the timescale at which decoherence saturates. 
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FIGURES 



f„(P) 




FIG. 1. A sketch of the mode function f n corresponding to a low-frequency mode of the elec- 
tromagnetic field, in the presence of a round superconducting wire (of radius R\). The sketch is 
meant to reflect two features: the main support of /„ is outside the superconductor, and the value 
of f n on the surface is much smaller than its typical value outside. 
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